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ABSTRACT 


Linear sound sweeps from a towed J15-3 sound source were collected at a moored 
VLA hydrophone array in the South China Sea during the ASIAEX experiment in May 
2001. Measured signals were filtered and pulse compressed. The processed data showed 
a high signal to noise ratio. Given an a priori chirp sonar survey, a two layer bottom 
“first guess” model was constructed. A broadband coupled-mode model was used to 
perform an exhaustive frequency variant sensitivity study of VLA pressures to changes in 
bottom properties as a basis for the geoacoustic inverse problem. Study results provided 
information on the observability of the various geoacoustic parameters and a procedure 
for the inversion. Matched field processing of the VLA data, using the same coupled¬ 
mode model, was then performed to calculate ambiguity diagrams from which 
geoacoustic parameter estimates were obtained. Since VLA pressure fields were not 
sensitive to changes in the sediment attenuation coefficient, a matched field technique 
that correlated the slope of modeled transmission loss to the negative slope of lOlog of 
the observed energy was performed in order to obtain estimates of the attenuation. These 
estimates showed a frequency dependent attenuation coefficient in the 50-600Hz 
frequency band. 
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I. 


INTRODUCTION 


A. ASIAEX OVERVIEW 

The Asian Seas International Acoustics Experiment (ASIAEX) was a joint 
scientific experiment between the United States, China, Korea and Taiwan (Chiu et al., 
2004). Institutions in the United States included the Naval Postgraduate School, Naval 
Research Laboratory and Woods Hole Oceanographic Institution. The study was multi¬ 
faceted but was primarily focused on exploring the predictability of low frequency sound 
propagation along/across the shelf break in the South China Sea and the effect of ocean 
variability upon that predictability (Ramp et al ., 2004). The prediction required a 
characterization of geoacoustic parameters and is the focus of this paper. 

The experiments were performed in April-August 2001 in the region of the South 
China Sea (Figure 1). Data used in this study was gathered on May 5 th 2001. This day 
was specifically chosen for its minimal water column variability (Chiu et al., 2004). 


South China Sea - ASIAEX 2001 



Figure 1. ASIAEX experiment location (shown in red box) relative to China 
and Taiwan coastlines situated on continental shelf/slope. Bathymetry is 

indicated by contours. 
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The acoustic source of this data was a J-15-3 broadband source towed behind the 
R/V Ocean Researcher III (OR3) at a speed of 2 knots and a depth of 3 meters. A variety 
of waveforms were broadcast including a linear frequency modulated (LFM) sweep 
ensemble used for this experiment (Figure 2). The LFM sweeps were broadcast 
continuously over 12 hours. Each thirty second ensemble was broken into three ten 
second sweeps with low, mid and high broadcasts of 50-200 Hz, 240-260 Hz, and 550- 
600 Hz hereafter referred to as LF, MF and HF respectively. A 10% cosine window was 
applied to the LFM waveform to properly ramp up amplitude. The frequencies chosen 
were specifically designed not to overlap in frequency with signals from moored sources 
that were broadcasting simultaneously (Turgut et al., 2003). 


5km LFM Specgram 



Figure 2. J-15 source signals used for geoacoustic inversion. Black arrows 

indicate low, mid and high frequency sweeps at 50-200 Hz, 240-260 Hz and 

550-600 Hz respectively. 

Received acoustic data was recorded on a moored hydrophone vertical line array 
(VLA) housed on an instrument sled in 124 meters of water. Hydrophone spacing on the 
VLA was set at 3.75 meters in order to adequately cover most of the water column with 
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its 16 elements. Data was continuously recorded over the entire duration of the towed- 
source operation with the ship track shown in Figure 3. 




longitude (deg) 


Figure 3. J-15 ship track (bottom) and range from VLA (top). Thin black 

arrows show the direction of ship movement. Green dot is the range (5km) 
where vertical measured data (black dot) for the matched field processing 
was gathered. Black x’s are the ranges used for horizontal data applied to 

the attenuation inversion only. 


B. BOTTOM MODEL 

The environmental bottom model structure used for both the sensitivity study and 
matched field processing was obtained by examining the chirp sonar survey imagery 
obtained by Steve Schock during ASIAEX (Schock, 2004). Visual analysis shows a 
simple two layer bottom structure, a thin sediment layer overlaying a thicker bottom 
layer, on the continental shelf and slope. This paper is focused only on the shelf as the 
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acoustic path from the towed source to the VLA never crossed more than a 20 meter 
change in bathymetry. 

While the chirp data provided the two layer structure of the bottom model (Figure 
4), core samples performed during ASIAEX provided a “first guess” sediment 
classification of the top meter of the bottom. Empirical and regressive equations obtained 
by Hamilton (1980) based on historical, world wide data were used to provide rough 
sound speed, density (Hamilton, 1977) and attenuation values given sediment type. 
Results from a previous inversion (Lin et al. 2001) in the same area but using a lower 
frequency broad band source from dynamite fishing were checked against the Hamilton 
data and, if reasonable, were used as first guess nominal values. 



Figure 4. Across shelf chirp sonar image of bottom bathymetry for ASIAEX 
location. The VLA receiver is located at the 124m isobath and is marked by 

a black pentagon. 


Bottom parameters in question, in each layer, are sound speed, attenuation and 
density. Given the two layer nature displayed by the chirp sonar, thin layer thickness 
becomes the seventh and final parameter, four in the thin layer and three in the bottom. 

C. THESIS OBJECTIVES 

The sonar operator with a detailed understanding of acoustic propagation and the 
environmental parameters that affect it carries a significant advantage over an adversary 
with all other factors equal. Therefore, the overall objective of this thesis is to provide 
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the warfighter with information critical to undersea warfare. As the typical naval 
operational environment migrates from open ocean “blue water” to littoral shallow water 
locales, an adequate characterization of the complex bottom environment is critical to the 
accurate prediction of acoustic propagation. The quantification of bottom characteristics 
in the experiment region via an inverse method is a means to help achieve the goal of 
warfighter support, in this case, to process the VLA data to obtain reliable estimates of 
the geoacoustic parameters. To ensure a sound scientific answer for the inversion, a 
sensitivity study is needed to determine the observability of acoustically significant 
parameters in the bottom, the quality of the inversion and the formulation for the inverse 
procedure. 

D. THESIS OUTLINE 

To accomplish the thesis objectives, the investigation into sensitivity of VLA 
pressure fields to bottom parameters was modeled first. Model runs were executed for 
each of the parameters in question. Results were cross-correlated to show vertical sound 
field sensitivity to each of the parameters, inversion quality and the information to 
develop a procedure for the inversion. 

With the sensitivity study completed, the model was again run to produce 
multiple VLA fields for matched field processing using the study findings as guidelines. 
Modeling was performed by the same coupled modes model. 

Details and results of the sensitivity study are presented in Chapter II. Chapter III 
explains the inverse technique and presents results of the geoacoustic parameter 
estimates. Conclusions and recommendations for further study are presented in the final 
chapter. 
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II. SENSITIVITY STUDY 


A. SENSITIVITY STUDY DETAILS 

The technique of matched field processing is not new to the scientific community 
or acousticians in general. For several decades this technique has been used extensively 
on a variety of candidate parameters for different regions (Lynch et al., 1990). What 
many papers fail to explore prior to performing matched field processing is the sensitivity 
of the modeled data output to changes in the candidate parameters. This is explored by 
the creation of a sensitivity diagram which shows the correlation coefficients of a large 
number of pressure fields generated by incremental changes in a candidate parameter. If 
little change is observed in output correlation coefficients given a large range of inputs, 
then the inversion is useless as a reasonable ambiguity diagram cannot be determined. 

Sensitivity of the model output also provides insight into the quality of the 
inversion. The narrowness of the peak in the sensitivity diagram describes the 
resolvability of the given candidate parameter and the accuracy that a matched field 
estimation can produce. The most important use for the sensitivity study is the 
information it provides allows for the development of an efficient inversion. Results 
from the sensitivity help determine which parameter is inverted for first, second and so 
on. 

The sensitivity study was performed using a broad band coupled modes model 
developed by Chiu (1996). This is the same model that was used for the matched field 
processing of the real data with the corresponding results presented in Chapter III. 

Seven parameters provided the scope of the sensitivity study. 

1. Thin layer sound velocity, c (mis) 

2. Bottom layer sound velocity, c (m/s) 

3. Thin layer density, p (kg / m 3 ) 

4. Bottom layer density, p (kg / m 3 ) 

5. Thin layer attenuation, a(dB/km /kHz) 
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6. Bottom layer attenuation a (dB/km/kHz) 


7. Thin layer thickness, h (m) 


The ranges and step increments for the candidate parameters were selected based 
on empirical results from Hamilton (1980) and a similar inversion performed in the same 
region but at lower frequency by Lin et al. (2001). The prior inversion was used as a 
basis for the parameter center values. Uncertainty ranges for each parameter were also 
garnered from the prior inversion but were expanded to ensure that all possible values 
were explored in the study. The empirical results obtained by Hamilton were used as a 
quality control. 

Each parameter was studied by holding the others constant while the one in 
question was varied over a given step increment. Bathymetry for the model was input as 
a flat bottom. Since the slope of the bathymetry was so gradual, a flat bottom simplified 
the model runs without affecting the results. Hydrophone receiver depths and source 
depth inputs matched the reported depths exactly. 

Frequency was varied for each parameter. Initially, the center frequency of each 
ASIAEX J-15 sweep was used (125, 250 and 575 Hz). However, the frequency range 
and increment changed for certain parameters as the initial center frequencies showed 
weak sensitivity or none at all. This was most noticeable in the MF and HF J-15 
frequency sweeps. Since the low sweep had the broadest frequency band and the greatest 
ability to penetrate the bottom, much information was gathered at these frequencies and 
necessitated a broader sensitivity study. 

Model output of the VFA pressure field was used for the sensitivity study. 
Pressure output was demeaned vertically and normalized to have the sum of the squares 
equal to unity. The output was then correlated to the conjugate transpose of pressure 
fields produced by incrementally changing the parameter to be studied. Calculated cross 
correlation coefficients were then contoured to give the sensitivity diagram. 
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B. STUDY RESULTS 

1. Thin Layer 

The thin layer is characterized by four parameters. It would be ideal if a 
frequency band existed that only showed sensitivity to the upper thin layer and none in 
the lower layer. If no frequency of that type could be found, the inversion would be 
much more computationally demanding. 

The first thin layer parameter examined was sound speed. Values varied from 
1550 to 1650 m/s with an increment of 5 m/s. Results showed weak sensitivity in the 
lower bounds of the low frequency sweep (50-75 Hz) as evidenced by the high 
correlation coefficient values with large standard deviations (wide peaks). Good 
sensitivity, shown by high correlation coefficients with small standard deviations, that 
improved with increasing frequency were observed for the rest of the frequency range 
(Figure 5). 

Values for density varied from 1520 to 1720 by 5 kg I in increments. Results 
were poor with no sensitivity showing for the thin layer (Figure 6). 

Attenuation values were varied between .1 and .35 by .025 dB/km/kHz 
increments. In prior inversions (Rajan et al., 1987) using a normal modes model, 
attenuation was not invertible so it was expected that attenuation fluctuation over a large 
range of inputs would not produce any sensitivity in the pressure fields. The results 
matched expectations (Figure 7). Attenuation would therefore need to be estimated via a 
different method. Instead of using VLA pressures, pressure vs. range outputs from the 
BBCM model were used to calculate TL. Figure 8 shows that slopes of the TL are 
sensitive too'. Vertical information is therefore not useful, but a range aperture is 
required for estimation of a . 

The final parameter for the thin layer was thickness. The input model values 
ranged from 2 to 20 meters every one meter. Results (Figure 9) showed moderate but 
acceptable sensitivity in the lower portions of the low sweep (50-75 Hz). Excellent 
sensitivity at the center frequency of the low sweep was observed. The mid to high 
frequency sweeps had poor sensitivity. 
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Sound velocity (m/s) Sound velocity (m/s) Sound velocity (m/s) Sound velocity (m/s) 


50Hz Thin Layer Sound Velocity Sensitivity 


75Hz Thin Layer Sound Velocity Sensitivity 





250Hz Thin Layer Sound Velocity Sensitivity 


575Hz Thin Layer Sound Velocity Sensitivity 


200Hz Thin Layer Sound Velocity Sensitivity 


125Hz Thin Layer Sound Velocity Sensitivity 


150Hz Thin Layer Sound Velocity Sensitivity 


175Hz Thin Layer Sound Velocity Sensitivity 



Sound velocity (m/s) 


Sound velocity (m/s) 





Figure 5. Thin Layer Sound Speed Sensitivity. Frequency increases from 50Hz 
at the top left to 575Hz at the bottom right. Frequency on diagrams increase 
left to right and top to bottom. Colorbar ranges from zero to one and 
represents correlation values. Sound speed increments range from 1550 to 

1650 m/s. 
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Density (gm/m 3 ) Density (gm/m 3 ) 


50Hz Thin Layer Density Sensitivity 75Hz Thin Layer Density Sensitivity 



1520 1540 1560 1580 1600 1620 1640 1660 1680 1700 1720 1520 1540 1560 1580 1600 1620 1640 1660 1680 1700 1720 

125Hz Thin Layer Density Sensitivity 250Hz Thin Layer Density Sensitivity 



1520 

1520 1540 1560 1580 1600 1620 1640 1660 1680 1700 1720 


1520 " 

1520 1540 1560 1580 1600 1620 1640 1660 1680 1700 1720 


575Hz Thin Layer Density Sensitivity 




Figure 6. Thin Layer Density Sensitivity. Frequency increases from 50Hz at the 
top left to 575Hz at the bottom. Frequency on diagrams increase left to right 
and top to bottom. Colorbar ranges from zero to one and represents 
correlation values. Lack of contours shows insensitivity of VLA pressure to 

bottom density fluctuations. 
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Attenuation (dB/km/kHz) Attenuation (dB/km/kHz) 


50Hz Thin Layer Attenuation Sensitivity 125Hz Thin Layer Attenuation Sensitivity 



Attenuation (dB/km/kHz) 

Figure 7. Thin Layer Attenuation Sensitivity. Frequency increases from 50Hz 
at the top left to 575Hz at the bottom. Frequency on diagrams increase left 
to right and top to bottom. Colorbar ranges from zero to one and represents 
correlation values. Lack of contours shows insensitivity of VLA pressure to 

bottom attenuation fluctuations 
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Transmission Loss Slopes for Modeled Data 125 Hz 



Transmission Loss Slopes for Modeled Data 250 Hz 



Transmission Loss Slopes for Modeled Data 575 Hz 



Source Range (km) with respect to 4km 


Figure 8. Thin layer attenuation sensitivity via TL slopes of modeled data. 
Subplots are for low, mid and high LFM sweeps. Frequencies increase from 
top (125Hz) to bottom (575Hz). Changing slopes show that TL slope is 
sensitive to a and that range aperture is required to estimate a . 
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Thickness (m) Thickness (m) 


50Hz Thin Layer Sensitivity 125Hz Thin Layer Sensitivity 



4 6 8 10 12 14 16 18 20 4 6 8 10 12 14 16 18 20 

75Hz Thin Layer Sensitivity 250Hz Thin Layer Sensitivity 
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575Hz Thin Layer Sensitivity 



4 6 8 10 12 14 16 18 20 

Thickness (m) 


Figure 9. Thin Layer Thickness Sensitivity. Frequency increases from 50Hz at 
the top left to 575Hz at the bottom. Frequency on diagrams increase left to 
right and top to bottom. Colorbar ranges from zero to one and represents 
correlation values. Thickness increments range from 2 to 20 meters. 
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2. Bottom Layer 

Visual inspection of the chirp sonar image showed the thin sediment layer 
covering the bottom layer to be at least 10 meters thick. This thickness is important 
number as penetration depth of the sound is a function of frequency. A 50 Hz signal (~30 
meter wavelength) can penetrate approximately 50 meters (Lin et al. 2001) with less 
penetration as frequency increases. The hypothesis was that the bottom layer would 
therefore be sensitive to the lower frequencies only. 

Parameter values were varied between 1610 and 1750 m/s by an increment of 5 
m/s. As expected, increased sensitivity was observed at the lower frequencies, especially 
in the lower portion of the low sweep (Figure 10). Poor or no sensitivity was seen above 
140 Hz (Figure 11) as the wavelengths became smaller than 11m. 

Attenuation and density showed similar poor sensitivity results in the bottom 
layer as in the thin layer. The graphical results were omitted in the interest of space and 
repetition for this paper. The results are identical to the thin layer density and attenuation 
graphics in Figures 6 and 8. 
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Figure 10. Bottom Layer Sound Speed Sensitivity (Low Frequencies). Frequency 
increases from 50Hz at the top left to 140Hz at the bottom right. Frequency 
on diagrams increase left to right and top to bottom. Colorbar ranges from 
zero to one and represents correlation values. Sound speed increments 

range from 1610 to 1750 m/s. 
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175Hz Thick Layer Sound Velocity Sensitivity 
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Figure 11. Bottom Layer Sound Speed Sensitivity (Low Frequencies). Frequency 
increases from 150Hz at the top left to 575Hz at the bottom. Frequency on 
diagrams increase left to right and top to bottom. Colorbar ranges from zero 
to one and represents correlation values. Sound speed increments range 

from 1610 to 1750 m/s. 
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III. DATA INVERSION 


A. MATCHED FIELD PROCESSING 

Matched field processing (MFP) has quickly become one of the most popular and 
widely used signal processing techniques for acoustics. The popularity is due to its 
ability to handle a wide range of problems that often confront acousticians in their 
investigations. Matched field is ideal for studies to determine parameters that affect the 
propagation of sound (Tolstoy, 1993), which is the goal of this paper. MFP also works 
well for localization type problems as well as model studies, again touched upon in this 
paper during the sensitivity study in Chapter II. 

Matched field does have its limitations. Noise is the enemy of MFP as is an 
ocean environment that is temporally variable during the period of data gathering 
(Tolstoy, 1993). . The J-15 propagation experiment occurred over 3 days, however, only 
the May 5 th data was used in this study as it was the most stable. 

As the field of MFP has progressed, numerous modifications have been added to 
improve results based on the variables in a given experiment. Processors have certain 
advantages over others. Resolution, sidelobes, and model errors are all handled 
differently with tradeoffs from processor to processor (Tolstoy, 1993). 

This study uses a linear (Bartlett) processor. It is the simplest and most 
understood of the processors available. It handles model inaccuracies and perturbations 
well and was a good choice for the study. Although sidelobes (regions of erroneous high 
processor power) are often associated with this processor (Tolstoy, 1993), they did not 
prove to be a problem. 

The concept behind the processor is simple. It correlates directly measured data 
with modeled data, in this case measured and modeled pressure fields. Processor power 
is defined as: 
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/ \ M / o o \ M 

C(x) = p (p p )p , 

where p M is the modeled VLA pressure field, p" is the observed VLA pressure 
field and p + is the conjugate transpose of p . Note that these are column vectors of 16x1 
elements and correspond to the 16 hydrophones in the VLA. 

Prior to entering the data into this linear correlator, both modeled and observed 
data was respectively normalized to have unit energy (sum of squares equal to one), 
except for the normalization for the attenuation discussed in detail below. 

Raw observed data was sampled every kilometer from 4 to 15 km. At each 
distance, 150 seconds of data were gathered (Figure 3.). As the LFM sweep ensemble 
was thirty seconds in duration, each range sample point had five LFM sweeps included. 
LF, MF and HF sweeps, each 10 seconds long, were combined into usable vectors. Data 
was then band-pass filtered and pulse compressed to further increase signal to noise ratio 
(SNR). Pulse compression used a synthetic replica signal to match each J-15 
transmission which included the 10% cosine taper used in the experiment. SNR was 
more than adequate after pulse compression for most frequencies. However, the cosine 
taper reduced the initial and final 5% of the observed signal which resulted in some SNR 
issues in the bottom layer inversion. The issue was resolved by choosing frequencies 
away from the edges of the linear sweep. 

B. ESTIMATOR DESIGN 

The final results from the sensitivity study showed that changes in three of the 
seven candidate parameters produce significant changes in the VLA sound pressures. 
The study also showed that relatively small (5%) changes in those parameters produced 
noticeable changes in pressure and would therefore produce an acceptable quality 
inversion. Most importantly, the study revealed information critical for formulation of 
the inverse procedure and is as follows: 


1. Invert for thin layer sound speed and thin layer thickness first. 
Frequencies used for the inversion should be between 140-200 Hz as 
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frequencies lower than that showed sensitivity to the bottom layer sound 
speed as well. 

2. Next, invert for bottom layer sound speed. Bandwidth acceptable for 
inversion is between 65-140 Hz. Signals with frequencies above 130 Hz 
have wavelengths too short to penetrate deep enough to feel the bottom 
layer changes. 

3. Density is not resolvable but can be determined using regression formulas 
by Hamilton (1980). 

4. Invert for attenuation using the modeled TL slopes. Center frequencies of 
each LFM sweep are to be used for the inversions. Vary attenuation 
values to obtain TL slopes and perform matched slope processing by 
comparing modeled TL slopes to the negative slope of lOlog of measured 
energy. 

C. THIN LAYER INVERSION 

The results from the sensitivity study proved that perturbations in two of the four 
candidate parameters, thin layer thickness and sound speed, resulted in appreciable 
changes in VLA pressure field differences. Therefore, only those two parameters could 
be inverted for using the acoustic pressure field as the base correlate input. 

Since sound speed is a parameter in both layers of the bottom, it was necessary to 
find a frequency that was sensitive to sounds speed in the upper layer but not in the 
bottom. The lowest frequencies (50-125 Hz) showed sensitivity to sound speed in both 
layers and therefore could not be used. The mid/low frequencies (140-200 Hz) showed 
good sensitivity to thin layer sound speed, no sensitivity to the bottom sound speed and 
were sensitive to thin layer thickness. Because of this, it was necessary to process both 
parameters at the same time. This was the preferred bandwidth used for the thin layer 
inversion. 

The modeled VLA pressure fields were created in a similar method used in the 
sensitivity study. However, for the inversion, actual bathymetry was used instead of a 
flat bottom (Wei, 2000). Since the inversion for the thin layer was attempting to discover 
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two parameters the model was executed a large number of times. Each thickness, from 2 
to 20 meters, in one meter increments, varied over the full range of candidate sound 
speeds, from 1550 to 1650 m/s by 5 m/s increments. This procedure was performed for 
each frequency in the preferred bandwidth in 5 Hz intervals. 

Observed data for the inversion was prepared in the manner described above for 
each discrete frequency to compare to its modeled counterpart. Although the movement 
of the ship in the time that it took for the J-15 to transmit from 50 to 200 Hz was 
minimal, and the distance traveled in a 5 Hz portion of that sweep was even smaller, it 
was critical to ensure accuracy of the field to be correlated. Any error in distance more 
than one wavelength produced significant degradation in the ambiguity surface. Results 
were very sensitive to distance. 

Ambiguity surfaces (Figure 12) were produced according to the above matched 
field methodology using a linear processor. Each surface was then analyzed to determine 
the maximum correlate value and its corresponding sound speed and thin layer thickness. 
Surfaces with a maximum correlate of less than .35 were discarded as inaccurate. Those 
meeting the correlate criteria threshold were then averaged using the following weighted 
mean formula. 


^c,hy=Yj 


w(c, h ) 


I 


w 


where c and h are the best sound speed and thin layer thickness estimates from the 
ambiguity surface ensemble and w is the corresponding correlate peak value for each 
frequency meeting the .35 correlation criteria. Variance was calculated from the 
weighted mean value. Results showed the thin layer to have a sound speed of 1575.68 
m/s and a standard deviation of 24.3 m/s. The thickness of the thin layer was 17.21 
meters with a standard deviation of 3.26 meters. 


It is worthy to note that ambiguity surfaces were produced for both coherent and 
non-coherent pressure fields. This was true for the sensitivity studies as well as the thin 
layer and bottom inversions. All graphics depicted in this paper are coherent. Surfaces 
constructed from magnitude only fields were too ambiguous to determine any of the 
parameters in question. 
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Figure 12. Thin layer ambiguity diagrams for both sound speed and layer 
thickness. Frequency increases from 140Hz at the top left to 180Hz at the 
bottom right. Frequency on diagrams increase left to right and top to 
bottom. Colorbars are not constant for each subplot with darkest red 
representing max correlation coefficient values for each diagram. Sound 
speed varies from 1550 to 1650 m/s while layer thickness ranges from 2 to 20 

meters. 
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C. BOTTOM LAYER INVERSION 

As the bottom layer only had one invertible parameter, the complexity of the 
model data matrix was much simpler than the thin layer. The matrix was constructed 
with model runs between 1590 and 1730 m/s with 5 m/s increments. Observed data was 
built in the same manner as the thin layer inversion. 

The frequency band used for the bottom inversion capped at 140 Hz based on lack 
of any sensitivity above that level. Frequencies below that level had varying degrees of 
sensitivity from 50-140Hz despite a clear diminishing sensitivity trend towards the high 
end of the bandwidth. 

Since the model runs for this inversion took exponentially less time to complete, 
more runs could be completed inside the frequency band used. Model runs were 
performed at 65, 79, 81, 86, 92, 102, 105, 112, 129, 137, and 141 Hz. The larger number 
of runs allowed for greater statistical significance of results. The irregular frequency 
interval is due to a lack of observed signal to noise ratio in the observed data. When 
regular spacing was used, ambiguity diagram correlate values were often very low due to 
a low SNR at some frequencies. In order to rectify the SNR problem, individual data 
points were chosen visually from FFT spectral values of pulse compressed raw data 
(Figure 13). 

Peak correlation values were chosen in the same method as the thin layer 
inversion except the threshold was .5 or higher. The same weighted sum was also 
employed. Sound speed in the lower layer was determined to be 1612.1 m/s with a sigma 
of 15.9 m/s (Figure 14). 
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Figure 13. Low frequency LFM amplitude spectrum after pulse compression 
showing obvious peaks in SNR that were used for bottom sound speed 

inversion in Figure 13. 
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Figure 14. Bottom layer sound speed frequency dependent ambiguity diagram 
with sound speed ranging from 1590 to 1730 m/s. Highest correlation 
coefficients are depicted in reds and oranges and correspond to estimated 

sound speeds on y axis. 
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D. ATTENUATION INVERSION 

As noted in the sensitivity study, VLA pressure fields did not respond to an 
incremented increase or decrease of attenuation values and were therefore impossible to 
invert for using only vertical data. 

The method that is valid and used here compared the slopes of the modeled depth 
integrated transmission loss (Figure 9) to the negative slope of lOlog of depth integrated 
measured energy (Figure 15) gathered from observed data. Modeled absolute 
transmission loss could not be compared directly to the observational data as the source 
level from the experiment J-15 sound source was not known. The two slopes can be 
correlated via matched field processing after some minor manipulation. Transmission 
loss is defined as: 

7X = -101og —, 

E 0 

Measured energy is defined as: 

L e =10 log E 

where E is the depth integrated energy and E 0 is the reference energy. The slopes 
of the model TL and —L e can be properly correlated to obtain high quality estimates of a . 

The modeled data vector ( p M ), which in this case consisted of multiple relative 
TL runs corresponding to increasing a values, was formed from the output of a 
broadband TL calculator developed by Chiu (1996). The data inputs for the TL 
calculator used pressure vs. range outputs from the same broadband coupled-mode model 
used in the other inversions for consistency. It is noteworthy to mention that the 
normalization procedure must use the larger of the modeled and observed sum of squares 
to make sure this is a matched slope processor. The same linear MLP as before was used 
except p° was replaced with —L e and p M was replaced with the modeled TL. 
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Peak correlation values pulled from the ambiguity diagrams (Figure 16.) produced 
by the matched field processor showed a non linear, frequency dependant attenuation 
(Figure 17). 


ASIAEx J15 LFM transmissions 



Figure 15. LFM measured negative energy (-lOlogE) slopes for low (50-200Hz), 
mid (240-260 Hz) and high (550-600 Hz) sweeps relative to energy at 4km. 


27 








Figure 16. 


125Hz Attenuation Ambiguity Surface 



250Hz Attenuation Ambiguity Surface 



575Hz Attenuation Ambiguity Surface 



Attenuation ambiguity diagrams. Top subplot is for 125Hz, middle is 
250Hz and bottom is 575Hz. X axis represents the attenuation value 
estimated from its corresponding TL slope. The y axis is the correlation 

coefficient value. 
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Figure 17. Attenuation rate versus frequency. Data set fit with a cubic spline 

curve. 
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IV. CONCLUSIONS 


Data collected during ASIAEX and the corresponding sensitivity study produced 
excellent results and provided insight into techniques for future geoacoustic inversions. 
Results show that matched field processing alone is not adequate to properly characterize 
bottom parameters. Investigations into sensitivity of the parameters must accompany 
matched field processing prior to the inversion to understand resolvability and accuracy. 

The sensitivity study produced one of the most important results of this thesis; the 
procedure describing the order to conduct the inversion is critical to obtaining efficient 
results. The sensitivity study also revealed the resolvability of the attenuation coefficient. 
Vertical pressure data is of no help to processing for attenuation. Horizontal data is 
useful and it can be estimated using matched field processing by matching the slopes of 
transmission loss to the negative slope of lOlog of observed energy. Preliminary results 
also suggest that attenuation rate is not linear with frequency as previously thought 
(Hamilton, 2001). That study suggested that attenuation rate = af 1 , regardless of 
sediment type. More research on this is recommended. 

The chirp sonar data was the most compelling data source to determine the 
layering of the bottom model. It clearly showed two layers of sediment with a thin layer 
overlying a thicker layer. Due to the layering of the bottom, and the thickness of the thin 
layer approaching 20 meters, the bottom layer becomes acoustically unimportant at 
higher frequencies. Penetration is not deep enough to discern the bottom layer 
parameters. 

This study focused on only one range (5km) from the experiment in order to 
invert for all of the parameters excluding attenuation which utilized horizontal received 
energy. Many more ranges could be explored as well as the all encompassing range 
dependent problem. 

As the overall goal of this thesis was Naval warfighter support, one has to ask 

how a study of this nature is processed to reach the operational models of a sonar 

operator on a surface ship or submarine. The gap between the operational Navy and 

research community needs to be tightened. Onboard acoustic models in use today have 
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little valid bottom input for shallow water available to the operator. Studies like this one 
have the potential to aide in the improvement of future operational models. 
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